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Abstract 

We consider the non-forward amphtude within the Heavy Quark Effective 
Theory. We show that one can obtain new information on the subleading cor- 
rections in l/rriQ. We illustrate the method by deriving new simple relations 
between the functions (,3(10) and A^^^w) and the sums r (w) 

(j = ^, |), that involve leading quantities, namely the Isgur-Wise functions 
T^^\w) and the level spacings A£'j"\ Our results follow because the non- 
forward amplitude B{vi) — > D^'^\v') B{vf) depends on three variables 
{wi,Wf,Wif) = {vi ■ v',Vf ■ v',Vi • Vf) independent in a certain domain, and 
we consider the zero recoil frontier {wjljw) where only a finite number of 

states contribute (^^) i^)' '^'^^^^ rules reduce to known results at 
w = 1, for A obtained by Voloshin, and for 1^3(1) obtained by Le Yaouanc 
et al. and by Uraltsev, and generalizes them to all values of w. We dis- 
cuss phenomenological applications of these results, in particular the check of 
Bakamjian-Thomas quarks models and the comparison with the QCD Sum 
Rules approach. 
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1 Introduction. 



In the leading order of the heavy quark expansion of QCD, Bjorken sum rule 
(SR) [11121 relates the slope of the elastic Isgur-Wise (IW) function ^{w), to the 
IW functions of the transition between the ground state = | and the = 
1^, 1^ excited states, Tyl{w), Tyl{w) at zero recoil w = 1 (n is a radial quantum 
number). This SR leads to the lower bound — 'C'(l) = p"^ ^ \- A new SR was 
formulated by Uraltsev in the heavy quark limit 3J, involving also Tyl{w), t!^^1{w), 
that implies, combined with Bjorken SR, the much stronger lower bound > |. 
A basic ingredient in deriving this bound is the consideration of the non-forward 
amplitude B(vi) D^'^\v') — > B{vf), allowing for general Vi, vj, v' and where B is 
a ground state meson. In refs. jH El we have developed, in the heavy quark limit 
of QCD, a manifestly covariant formalism within the Operator Product Expansion 
(OPE), using the matrix representation for the whole tower of heavy meson states 
[71. We did recover Uraltsev SR plus a general class of SR that allow to bound also 
higher derivatives of the IW function. In particular, we found two bounds for the 
curvature ^"(1) = cr^ in terms of p^, that imply cr^ > jf- 

The object of the present paper is to extend the formalism to IW functions at 
subleading order in I/ttiq. 

The general SR obtained from the OPE can be written in the compact way [1] 

LnadronsiWi, Wf, Wif) = RoPE{Wi,Wf, Wif) (1) 

where the l.h.s. is the sum over the intermediate D states, while the r.h.s. is the 
OPE counterpart. Using the trace formalism [HI, this expression writes, in the heavy 
quark limit ^ : 



D=P,V n 



J2 E Tr [Bf{vf)TfD^^''\v')] Tr [D{v')T,B^{vi)] ^^"^K)^^''^^^/) + Excited states 

= -2^K/) Tr \Bf{vj)TfP[T,B,{v,)] (2) 



where 



and 



Wi = Vi-v' Wf = Vfv' Wif = Vi-Vf (3) 



n - (4) 
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is the positive energy projector on the intermediate c quark, and we assume that 
the IW functions are real. The ground state B meson can be either a pseudoscalar 
or a vector, the doublet | . The heavy quark currents considered in the previous 
expression are 

(5) 

B{v), D{v) are the 4x4 matrices representing the B, D states [HI, and B = 7^5+7° 
denotes the Dirac conjugate matrix. The domain for the variables {wi,Wf,Wif) is 

Wi > 1 Wf >1 

WiWf - ^{wf - l){w) - 1) < Wif < WiWf + ^{wf - l){w) - 1) . (6) 

We will now consider l/mg corrections to the heavy quark limit SR (0). 

The paper is organized as follows. In Section 2 we set the problem of obtaining 
sum rules involving subleading quantities in l/mq within the OPE. In Section 3 
we make explicit the formalism of the corrections in l/rrib, the 6-quark being the 
external quark, using the formalisms of Falk and Neubert ~9j and of Leibovich et 
al. [TU] to parametrize the | — > | and | 1^? §''' form factors, that we extend 
to the I ^1 transitions. In Sections 4 and 5 we write down the SR obtained 
respectively if the initial and final meson is a pseudoscalar i? or a vector B*. We 
factorize polynomials in the variables {wi, Wf,Wif) that allow to obtain simple results 
going to the interesting frontier {w, 1, w) of the domain ©. We generalize our results 
to any ^ transition. From the obtained SR we get enough information to 

write our fundamental results for the subleading quantities A^{w) and ^3(10) in 
terms of leading quantities in Section 6. In Section 7 we use as input the results 
of the Bakamjian-Thomas class of quark models - that satisfy all the necessary 
properties in the heavy quark limit - to obtain phenomenologically useful results, 
that appear to be consistent, in our quite different approach, with the QCD Sum 
Rules. Finally, in Section 8 we conclude and set up the program that remains to 
be pursued. In Appendix A we demonstrate the identity between two subleading 
parameters defined by Falk and Neubert ^9^ and by Uraltsev 3J. In Appendix B we 
discuss the experimental situation of the leading P-wave IW functions Ty2{w) and 
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2 The Operator Product Expansion and the 
corrections at first order in l/mg. 

Our starting point [TT] is the T-product 

TM =ijd^x e-'''-^ < 5(pj)|T[JXO)J,(x)]|5(pO > (7) 

where Jf{x), Ji{y) are the currents (the convenient notation for the subindices i, f 
will appear clear below) : 

Jf{x) = b{x)Tfc{x) My) = c{y)TMy) (8) 

and Pi is in general different from pf. 

Inserting in this expression hadronic intermediate states, a;*^ < receives con- 
tributions from the direct channel with hadrons with a single heavy quark c, while 
x° > receives contributions from hadrons with bcb quarks, the Z diagrams : 

X, q'' + hi - J^Xc + 

- E (2.r - Pf + Px,J ^ ^/l-^-'yi-^"^ >< J;(0)|B. > 

We will consider the following limit 

rric':^ rrib':^ Aqcd ■ (10) 

The difference between the two energy denominators is large 

qO-Ef + Ex-^,, - (g° + Ex) ~ 2m, . (11) 

Therefore, we can in this limit neglect the second term, and we will consider the 
imaginary part of the direct diagram, the first term in (jH)), the piece proportional to 

5((f + E,-Ex}, . (12) 

One can see this point otherwise. The two cuts corresponding to the two terms in 
(P) are widely separated, and one can isolate the imaginary part of the first term by 
a suitable integration contour in the g° complex plane. Notice that one can choose 
such that there is a left-hand cut, even in the conditions ()10|) . This means that 



is of the order of and — is fixed, of the order mi,. Our conditions are, in 
short, as follows : 

Aqcd <^ nib nic - q° <^ q° rric , (13) 

consistent with ()12j). To summarize, we are considering the heavy quark limit for 
the c quark, but we allow for a large finite mass for the b quark. 

Unlike the case of the forward amplitude {pi = Pf), the imaginary part of the 
direct diagram in Q will not be related to a positive definite absorptive part, because 
we are in the more general case of the non-forward amplitude. However, we are 
allowed to consider this imaginary part. 

In the conditions (|T!?jl . or choosing the suitable integration contour [T2j [IH!, we 
can write therefore, integrating over q^ 

7>t(q) = E(2^)' (q + p. - qxj < Bj\Jj{o)\x, > < x^iom > . (i4) 

Finally, integrating over qx^ and defining v' = one gets 

r;f = ^ < Bfivf)\JfiO)\D^iv') > < D^iv')\MO)mv,) > (15) 

where we have denoted by Dn{v') the charmed intermediate states. 

The T-product matrix element Tfi{q) (jZj) is given, alternatively, in terms of 
quarks and gluons, by the expression 

Tj,{q) = - j (fx e-"^-^ < B{pf)\mTjS,{Q,x)T,h{x)\B{p,) > (16) 

where Sc(0,x) is the c quark propagator in the background of the soft gluon field 

m- 

Since we are considering the absorptive part in the c heavy quark limit of the 
direct graph in 0, this quantity can be then identified with (fTB|) where Sc{x, 0) is 
replaced by the following expression J3] 

S,iO,x) ^ e*'"^^'-^ <l>AO,x]D,,ix) (17) 

where Dy/ (x) is the cut free propagator of a heavy quark 

-|--5(fc ■ vy^- = p; ^5\x - v't) (18) 
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with the positive energy projector defined by 



n = ^(i+^')- (19) 



The eikonal phase $t,'[0, x] in (fTTjl corresponds to the propagation of the c quark 
from the point x = v't to the point 0, that is given by 

[0, v't] = P exp (^-i ds v' ■ A{v's)^ . (20) 

This quantity takes care of the dynamics of the soft gluons in HQET along the 
classical path x = v't. 

Inserting (fT7j) - ()19|) into we obtain 

r dt 

Tf^-(q) = / d^x e-^(™.')- / !^s\x - v't) 



-oo 27r 

< B{pf)\b{0)TfP'^<^,>[0,x]Tib{x)\B{pi) > + 0(l/m,) . (21) 
Integrating over x in (j2T| and making explicit (j20|l . 

•^^ J-oo 2tt 

< 5(j9/)|6(0)r/P;Pexp {^-i f^ds v' ■ A{v's)j Vih{v't)\B{pi) > 

+ 0(l/me) . (22) 

Performing first the integration over q'^ one obtains simply 6{v'^t), that forces 
t = 0, and making the trivial integration over t one obtains finally the OPE matrix 
element that must be identified with (IT5ll : 



^/t = < Bipj)mTj^TMO)\Bip.) > + 0(l/mj . (23) 
Therefore, we end up with the sum rule 

y: < Bjivf)\jfmD4v') > < D„{v')\Mo)mik) > 

= < 5(^;^)|6(0)^^^^,6(0)|5(^;,) > + 0(l/m,) (24) 

that is valid for all powers of an expansion in l/nih, but only to leading order in 
l/rric. 
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On the other hand, making use of the HQET equations of motion, the field b{x) 
in ()24|1 can be decomposed into upper and lower components as follows ^H! 

b{x) = e-''""^-^ I 1 + =, 0] K{x) (25) 

V 2mb + iv D ) 

where the second term corresponds to the lower components and can be expanded 
in a series in powers of -D^/m^, and v is an arbitrary four- velocity. 

Taking into account the normalization of the states in the trace formalism and 
the sign convention for the matrix elements, we recover, in the heavy quark limit 
(neglecting the term in l/2mfe), the master formula (2) obtained in ref. |3]. 

Including the first order in 1/m;,, the sum rule reads 

E < B^{v^)\Jm^n{v') > < Dr.iv')\MQ)mv.) > 
1 ..^ ,„.r, .t..„ l+l l+l, 



<B{pf)\K^{0)[{-z^)Tf^T, + Tj^TSl/)) KXm{Pr)> 



+ 0{l/m,) + 0{l/ml) . (26) 
Therefore, in the OPE side we have, besides the leading dimension 3 operator 

O^'^ = K^TfP:^T,K^ (27) 

the dimension 4 operator 

0(4) = ^, [{-tp)TfP[T, + TfP'.TS^)] . (28) 

In the SR we have to compute the l.h.s. including terms of order l/2mf,. These 
terms have been parametrized by Falk and Neubert for the | doublet and by 
Leibovich et al. for the transitions between the ground state | and the |^ 
excited states. These l/m^ corrections are of two classes : perturbations of the 
current, and perturbations of the Lagrangian (kinetic and magnetic). This will be 
our guideline to compute the l.h.s. of the SR, although we will consider the whole 
tower of excited states. 

A remark is in order here, that was already made in ref. [T^. Had we taken 



higher moments of the form / dq^ {q^)"'T?^^ {q^) [n > 0), instead of the lowest one 
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n = 0, the integration over g° that leads to the simple sum rules ()24|) or ()26p would 
involve higher dimension operators, giving a whole tower of sum rules JSj, even 
in the leading heavy quark limit. Our point of view in this paper is different. We 
consider the lowest moment n = 0, while we expand in powers of l/mt, keeping the 
first order in this parameter. 

Concerning the OPE side in ()2fj|l . the dimension 4 operator O*^^^ ()28j) is nothing 
else but the l/m^ perturbation of the heavy current O^^^ = h^^T fP'jVih^^ since this 
operator, containing the Dirac matrix T fP'jTi between heavy quark fields, can be 
considered as a heavy quark current. Indeed, following Falk and Neubert, the 1/mi, 
perturbation of any heavy quark current hyfVhm is given by 



However, this perturbation of the current does not exhaust all perturbations in 
l/rrih. Indeed, we need also to compute the perturbation of the initial and final 
wave functions \Bi{vi) >, \Bf{vf) > due to the kinetic and magnetic perturbations 
of the Lagrangian. This can be done easily following also the prescriptions of Falk 
and Neubert to compute these corrections in l/rrib for the leading matrix element 
< Bf{vf)\hvfT fP[_Tihm\Bi{vi) >, as we will see below. 

3 Setting the formalism for the calculation of the 
corrections in l/mb- 

Considering B or B* initial and final mesons, we can perturb the SR ((2)) by 
l/rric and l/rrib terms. The perturbation of the r.h.s. is parametrized by six new 
subleading IW functions concerning the ground state | , denoted by Li{w) {i = 
1, ■ ■ ■ 6), in the notation of Falk and Neubert p. 

As for the l.h.s., considering for the moment as intermediate D states the mul- 
tiplets I , 1^, we have three types of matrix elements 




(29) 



<Dii-)iv')\cTb\Biv)> 
<Dii^){v')\cTb\B{v)> 
<£)(§+) {v')\cTb\B{v) > 



(30) 
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The corrections in l/rri}, or l/m^ to the first matrix element are given by the same 
ground state subleading IW functions Li{w) (i = 1, ■ ■ -6), while the 0{l/mb) and 
0{l/m^ corrections to the matrix elements B ^ D (|^); D (f"*^) have been care- 
fully studied by Leibovich, Ligeti, Steward and Wise [TU] , and result in a number of 
new subleading IW functions. All these corrections are of two types, perturbations 
of the heavy quark current, and perturbations of the Lagrangian. 

Moreover, since, as pointed out by Leibovich et al. [10, Section VI], the states 
D contribute also to zero recoil at order l/mg, we will consider the contri- 

bution of the matrix elements 

<D{l-,j-){v')\cVh\B{v)> (J =1,2). (31) 

We will show that these contributions do not spoil the simple result presented 
below, that can be expressed only in terms of the leading IW functions tii2{w) and 
T-i/2{w). We argue also that higher intermediate states do not contribute. 

Let us again underline that we will not take into account radiative hard gluon 
corrections, as computed in jT2] for Bjorken SR, in for Uraltsev SR and in jT3] 
for our SR concerning the curvature of the IW function [H]. 

We begin with the general SR in the heavy quark limit Q and perturb the 
heavy quark limit matrix elements with l/m^ and l/m^ corrections. The general 
expression could then be written, making explicit the leading and the l/m-c and 
l/mft parts : 

Go{wi,Wf, Wif) + Eo{wi, Wf, Wif) + [Gb{wi, Wf, Wif) + Eb{wi,Wf, Wif)] 

f*^"^*^^" ^ Ec{wi, Wf, Wif)] 

= Ro{wi, Wf, Wif) + -^Rb{wi, Wf, Wif) + -^Rc{wi, Wf, Wif) (32) 

where the subindex means the heavy quark limit, while the subindex 6 or c corre- 
spond to the subleading corrections in l/nif, or l/rric, and G or E mean, respectively, 
ground state or excited state contributions. 
In the heavy quark limit, one has 

Go{wi,Wf, Wif) + Eo{wi, Wf, Wif) = Ro{wi, Wf, Wif) (33) 

that leads to equation (j21) and to the results quoted above PQ-IH]- 



In expression ()32j) we can vary rrih and as independent parameters and obtain 
new SR for the subleading quantities. 

To obtain information on the l/m^ corrections, it is relatively simple to proceed 
as follows. We will assume the formal limit of Section 2 : 

mc > mfe > Aqcd (34) 

and perturb both sides of the SR ()33|) by 1/ rrify terms. This heuristic procedure gives 
the same results as the method demonstrated in Section 2. 

In this limit, since the parameter l/rrib can be varied at will, one obtains the 
relation 

Gb{wi, Wf, Wif) + Eb{wi,Wf, Wif) = Rb{wi, Wf, Wij) . (35) 

One can compute Gb{wi,Wf,Wif) and Eb{wi,Wf,Wif) using respectively the for- 
malism of Falk and Neubert jH] and the one of Leibovich et al. |10j, and obtain SR 
for the different subleading IW functions Li{w) (i = 1, ■ ■ -6). 

Of course, one can obtain SR by taking the opposite limit rrib ^ rric, that must be 
consistent with the preceding ones. In ref. ^T] we did adopt the Shifman-Voloshin 
limit [IH] rrib, rric^ nib — rric^ ^qcd for the forward amplitude. 



To be explicit, let us define these functions from the current matrix elements, 
following the notation of Falk and Neubert [Hj : 



< D{v')\Q'TQ\B{v) > = -^{w) Tr [d{v')TB{v) 
— Tr {d{v')T [P+L+{v, v') + P^L_{v, v')]} 



2mb 

in the formal limit rric ^ mb fl34|) that we adopt here. 

The 4x4 matrices write, respectively, for pseudoscalar and vector mesons : 



(36) 



M{v) = P4v) 4 

while the subleading l/rrib functions are for pseudoscalar and vector mesons : 



(37) 



P4v)L4v,v') + P4v)L4v,v') = [L,iw)P4v) + L,iw)P4v)]i-^5) (3J 
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P+{v) l^M^) + i^v ■ v') L,{w)] + P_{v) i^M^) + i^v ■ v') L,{w)] (39) 

where w = v ■ v' . 

The matrix elements to excited states write ^Ul 

< D (1+) {v')\cTh\B{v) > ^ V3 T^/2{w) Tr [v.D" {v')TB{v) 

+ ^{Tr [s';Id\v')V^^B{v)\ +r]t} Tr [vJD" {v')VB{v 



2mh 
+ Tr 



R^l,D\v')VP^{v)^a-'B{v)\ } 



where 



< D (1+) {v')\cVh\B{v) > ^ 2n/2iw) Tr D{v')TB{v 



+ 



+ Tr 



-{Tr [S^^^D{v')r^'B{v)] + x'ul Tr [d{v')TB{v) 
Rfp{v')VP+{v)ia^PB{v)\ } 



(40) 



D-,^iv') = P+iv')e7^, 
Di+(t;')=P+(t;')<'757. 

Do+iv') = P+iv') . (41) 

The notations 5*^^^ S^^^ denote the perturbations to the current, and rjl^^, xtl ^^"^ 
R^^ap and R^^p denote respectively the kinetic and the magnetic perturbations to the 
Lagrangian. In the preceding relations pn)) B{v) can be a pseudoscalar or a vector. 

Expanded in terms of Lorentz covariant factors and subleading IW functions, 
these tensor quantities read [121 • 



S. 



crA 



^ {w)vx + T^''^ {w)v'x + (w)7a + rf ^ {w)g^x 

^S = xS'^H7.7/3 + xf(^)<7/3. (42) 

The IW functions relevant to the current perturbation are not independent, due to 
the equations of motion : 



{b)i 



-ib), 
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.(b) 



(6)/ 



(43) 



and, at zero recoil, one has 



rf(l) 
Cf(l) 



V3 A^3/2 r3/2(l) 
-AEi/2 ri/2(l) 



(44) 



where a radial quantum number n is implicit and A.E1/2 are the mass differ- 

ences between the excited states and the ground state. 

Since, as pointed out above, we will also consider the intermediate states D (^| , , 
let us give the relevant formulae, parallel to (jlUll - PH) 

< D (§-) {v')\cTb\B{v) > = V3 (T3/2(«;) Tr [v^D'" {v')TB{v 
+ :7^{Tr [T!!'^D''iv')T^'Biv)] + p£ Tr [v^W {v')TBiv) 
+ Tr [v}%D''iv')TP+iv)ta''^Biv)] } 
Z}^(t;') = I^r+K)(-75) 



T, 



{w)vx + {w)v'^ + o-f ^ (w)7a + af^ {w)g„x 
af(l) = v^AE(r)a3/2(l). 



(45) 



At zero recoil, Luke's theorem ^H] imposes 

Li(l) = L2(1) = 



(46) 



while it can be shown that -^^4(1), -^^5(1), -^6(1) are not linearly independent jH], and 
are related to two quantities, namely 

A = niB — nib = niD — Trie (47) 
and the quantity called ^3(1) by Falk and Neubert or E by Uraltsev jH] : 

L4(l) = -A + 2e3(l) 
L,{1) = -A 

Lg(l) = -A-e3(l) • (48) 
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We demonstrate the identity ^3(1) = S in Appendix A. Considering the forward 
amphtude, i.e. taking Wij = 1, two SR can be obtained for subleading corrections 
at zero recoil, as we will see below : 

A = 2 E ^E[%\ri^l{l)\^ + 4 E AE^^l\4^l (1) ^ 

n 71 

6(l)^S = 25:AE(;J|r3(;)(l)p-2EA<^|r|;)(l)p (49) 

n n 

where Ai?]"^ and Tj"''(l) (j = |, |) are the corresponding level spacings and tran- 
sition IW functions between the ground state | and the P-wave states and 
I"*". The first SR is Voloshin SR ^j, and the second one was discovered by A. Le 
Yaouanc et al. |llj and by Uraltsev We have adopted the notation of Isgur and 
Wise for the transition IW functions j2j. 



4 B Meson sum rule. 



We take as initial and final states the ground state pseudoscalar meson at differ- 
ent four-velocities B{vi) and B{vf) and, as in the axial currents aligned along 
the corresponding four- velocities, Fj = jjf^^ and T/ = fib- Then, the subleading 
SR dSSI) writes : 



+ 



[WiWf 



if 



f-kwi-iM-i) 



3/2 



Wif)\/?) ayl{wi) Gy!2{'^f) + Higher states 



+ [WiWf 

= - [2Li{wif){l + Wif -Wi- Wf) - 2L4^{wif){l - Wif)] 



/)} 



(50) 



The first, second, third and fourth term in the l.h.s. of (j5(J|) correspond to the 
Q (2^' 2^) ^^'^ (1 5 1~) intermediate D states. No other states 

with J < f appear in the l.h.s. because of the number of 75 matrices involved in 



the traces over Dirac matrices. 

In equation (j^UI) we have made explicit the subleading ^ 

-in) 



elastic functions 



(j — l'^) and we have factorized, when possible, for the contributions | 



2 ' 2 



and I , polynomials in {wi,Wf,Wif) that vanish at the frontier (w,l,w) of 
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the domain The IW functions Tyl{w), Tyl{w) and cr^liw) do not vanish at 
w = 1. Similarly, the complicated functions F^'jl{w), Fyl{w) and Gy\{w) are given 
in terms of the form factors defined in Section 3 and do not vanish in general for 
w = 1. 

It is not necessary to give all the explicit expressions of these functions since we 
are interested in the frontier of the domain {wi,Wf,Wif) = {w, l,w) and, as we will 
see below, only F^jl{l) will contribute. 

There are two crucial features in expression (jSHJ. First, the appearance of the 
subleading functions Li{wif), L/^iwij) in the r.h.s., since we consider the whole al- 
lowed domain for the variables {wi, Wf, Wij). Second, the polynomials in {wi, Wf, Wif), 
that result from the sum over the spin J = 1, 2 polarizations 0] : 



4' 4' Vf^Vi^ =WiWf - Wif 

X 



^4^^*^^" 4^^^^" Vff^Vf^VipVi^ = {wiWf - Wiff - \{w'l - \){w) - 1) (51) 

where A runs over the 2 J + 1 polarizations. 

These polynomials will imply the vanishing of the corresponding contributions 
at {wi,Wf,Wif) = {w,l,w). This will occur also for higher intermediate states, 
because one obtains, in all generality |3], for the projector on the polarization tensor 
of a particle of integer spin J, contracted with Vi and Vf four- velocities : 

rpUl--Uj,tll --flj _ ^2^(X)*Ul—Uj ^(A)^ll---^Jj {52) 

A 

Vfiy-^ ■ ■ ■ Vfyj T^/ "'''^^ fj^j ■ ■ ■ f j^j = 

that vanishes for J > at {wi, Wf,Wif) = {w, l,w). 
Therefore, at the frontier 

(wi^Wf, Wif) (w, 1, w) (54) 
the SR (plj) will write, very simply, dividing by a factor {w — 1) 

Uw) = J:A/1{w)fIJI{1). (55) 

n 
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We only need the functions Fyl{w) for w = 1. The calculation gives, for all w, 



Fijliw) = (l-w) [Cf )(")(^) + (1 + 2w)C?^-\w) + + Qx^^'-^w) 

-2(1 + w)x$^^''\w)] + 2(1 + 2w)Cf ^"^(w) (56) 



and for w 



F(;)(i) = 6cf)(")(i) 



and from the relation ()44p we obtain finally 



'1/2' 



w) 



(57) 



(5J 



5 B* meson sum rule. 

We take as initial and final states the ground state vector meson at different 
four- velocities B*^^^\vi) and B*^^f\vf) and we adopt the particular case of the B* 
polarizations (see appendix A of ref. |l4|) : 



Si 



Vf - WifVi 



if 



1 



Vi - WifVf 



(59) 



that satisfy Ei ■ Vi = Ef ■ Vf = and ef = e^ 



-1. With the definitions ()59j) one has 



Ei ■ Ef = Wif, but we can change one global sign in to make Ei ■ Ef = —Wif and 
therefore Ei ■ Ef ^ —1 when Vi ^ Vf. The sum rules, being linear in Ei and in e^, do 
not depend on this overall sign. 

Then, performing the relevant traces, the subleading SR writes 

n 

-{l-Wf)Lt\wf) + {l+Wf)Lt\wf) 

+ 2T^;;l{Wi) [(e) ■ E,) K[''\wf,W,) + (e, ■ Vf) {e} ■ v') Kt\wf,w,) 

+ {e} ■ Vi) {e, ■ v') Kt\wf) + (e, ■ v') [e) ■ v') Kt\wf, w„ Wif) 

+ VS T3/2{Wi) (e} ■ Ei) S[''\wf,Wi) + (e) ■ v') {Ei ■ V f) Si''\wf,Wi) 
+ (e} ■ Vi) {Ei ■ v') S3'\wf,Wi) + {Ei ■ v') (e} ■ v') St\w f ,Wi,Wif) + {i 

+ ^til(wi)T'^'^\wf) [{w,f - w^Wf)Tr f f ^j:*^^ Tr fij^il^ _ 
+ Tr \j)fi)i)i^il^ Tr U^ij'f^f^^ 



/)} 
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+ V3al}l{wi)U'^''\wf,Wi) Tr [Y^'^^f^^,] Tr 

+ Vff,e*f^Tl^,''''"^VipVia{ei ■ v')V^''\wf, Wi) 

+ VfpVfuTl^,''''"'vip6ia{e*f ■ v')V:^''\wf,Wi) 

+ Vff,Vf^Tl^,''''"'vipVi^{e*f ■ v'){ei ■ v')v}''\wf,Wi) 

+ Higher intermediate states + ^ /)| 

= -{ - 2 (^e*f Si) {1-Wi-Wf + Wif) + (e) ■ v'^ {si ■ Vf) 

+ is) ■ {si ■ v') - (e*f ■ Vi'j {si ■ Vf) L2{wif) 

- {e, ■ Vf) (^e*f ■ v'^ {w,f - 1) + (e} ■ Vi^ {si ■ v') {w,f - 1) 

- {si ■ Vf) {e*f ■ v^ {wi + Wf -2) L-i{wif) 

+ 2 {e*f (1 -Wif) + [e) ■ v^ {si ■ Vf) L^{wif) 

+ [- {si ■ Vf) (e) ■ v'^ {w,f + 1) - (^e*f ■ w^) {e, ■ v') {w,f + 1) 

+ {si ■ Vf) {^e} ■ {wi + Wf-2) Leiwif)} 



(60) 



In the Lh.s. of equation ()60|) . the first term corresponds to the intermediate 
states Q , 0^^ + (^1 , 1^^ (both spins contribute due to the fact that one has more 
four- vectors and a 75 in the traces than in the pseudoscalar case), the second, third 



and fourth terms correspond to the contributions 



1+ 1+ 
2 ' 



' V2 ' 



and 



3 + 
2 ' ■ 



and the fifth and sixth terms correspond to the contributions (^| , and , 2^^ . 
We have made exphcit the subleading | | functions L^"^ (j = 2, 3, 5, 6) and 
we have kept the exphcit dependence on the initial and final polarizations Ei, Sf. 
This allows to factorize, when possible, for the contributions \ I"*") f'*' ^"^^ ^ 



2 ' 



polynomials in {wi,Wf,Wif) that vanish at the frontier {w, l,w). The IW functions 



Tyl{w), T3/2(^) ^^"^ '^3/2 (^) ^^^^ ^'^ vanish in general for w 
complicated functions /sfj"^ (j = l,---4), 5*]"^ (j = l,---4), T^"), f/(") and Vj 
{j = 1, ■ ■ - 4) can be computed in terms of the form factors defined in Section 3. 

The tensor T^^,"''"^ , that appears in the l.h.s. of eq. dEDl), is given in terms of the 
J = 2 polarization tensors by the expression 



1, and the 



rin) 



rpflU, 



pa 



(\)^u {X)po- 



(61) 



In order to see clearly which terms survive at the frontier {wi, Wf, Wif) = {w, 1, w), 
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it is not necessary to go to the details that we have given in Section 4 for the sake 
of clarity. It is enough to realize that the limit 

{wi,Wf,Wif) ^ {w,l,w) (62) 

corresponds to the limit 

f/ — > f ' , Vi ^ V (63) 

{v ■ v' = w) and make use of the orthogonality conditions between the interme- 
diate states polarization tensors and v'. Explicit calculations confirm this simple 
argument. 

In the limit dHBI), we have 



EfV 



Si ■ V 



Sj ■ Vi 



w 




-V w'^ — 1 



Si-Vf- 



1 



Tr 
Tr 

T'^'^'P"' „^ r* 7,. C-. 







rpfJ. 



fiu,pcr 







. 



VfpVfiyVipEia 



VfpVfyVipVi^ 



(64) 



The last limits follow from the orthogonality condition T^^'P"v'p = T^^'P"v'^ = 0. 
Notice that the "Higher intermediate states" contributions in equation ()60p. 
similarly to the four last expressions (j6H) and to the general expression (j53|l . due to 
the symmetry in (z/i, z/2, ■ ■ ■ z/j) and in (/xi, /i2, ■ ■ ■ fij) and the linearity in Ei and e*^, 
will be proportional to the following quantities 



rpUil>2-l'j,p.lP2-P.J „, „, 



rpUiU2-Uj,p.itJ.2-pj . 
rpUlV2---y,J,P.lP2---p.J 



rT.UlU2-Uj,fJ.ip2-pj 



' ^ fuj 

















(65) 
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The tensor t^^^'^^'"'^''''^^'^^'"'^"' is given, in terms of the polarization tensor of an inter- 
mediate state of spin J, by expression (jK^ . and the last quantity in (jHHjl is given by 
the polynomial in {wi,Wf,Wif) (jSSl)- 

By the same argument as before, due to the orthogonality conditions 



rpl'lV2---Uj,^ll^l2-^lJ I _ n 



(fc = 1 ■ ■ ■ J) 



(66) 



all the quantities (jHSl) go to in the limit Vf ^ v' 

Therefore, in the limit {wi,Wf, Wif) — » {w, 1, w), the SR taking into account 
its symmetry in z ^ / and the asymmetry in z ^ / of the limit becomes the 

much simpler expression : 

n 

+ 2r(;)(l) [wKi-\w, 1) + {w' - l)Ki-\w, 1)] } 



+ 



+ 4';liw) [wSt'\w, 1) + {w^ - l)St\w, 1)] } 
= -2{w - 1)L5H + 2{w^ - 1)Lq{w) . 



(67) 



Therefore, we have only to compute the functions Ki"'\wi,Wf), {wi,Wf), 
K^"'\wi), Si^\wi,Wf), S2"'\wi,Wf) and S^\wi,Wf). The exphcit calculation gives, 
for general {wi,Wf), 



Ki''\w,,Wf) = -(m, - l){wf - 1) Cr'^'iwi) + {2w, + i)cr^"^(^.) 



(b)(n). 



Win) I 



Kt\w,,wf) = -(1 - wf)[cr''>{w,) + {2w, + ixr^^iw.) - 2cr'"{w,) 

+xti''\w.)-2x?^''\w.) + 2x?^"\w,)_ 

Kt\wi) = -(1 - w,) [cf)("H«;.) + i2w, + l)C?^''\w,) 



Win) I 



.{b){n), 



(68) 



and 



1. 2 



in), 



St\w.,Wf) = -^{2-Wf)Gr{w.) 



in) I 



(69) 



where in the preceding equations 



G^^\w.) = (2 - ^.)rf + (2i^. + 1)(2 - ^.)rf ("^K) + 2ti;,rf ("V.) 

-2u;,rf (-)(^,) + (2 - - 2(2 - w,)vf^^-\w,) 

+4K - lW2^^''\w,) + 3{w, - 2)r/f ("^K) (70) 

and we get, for the quantities needed in equation ^7\i. 

ir{")(l,^) = -2(l-«;)Cp)(l) 



4 



-WW, 



5(")(l,w) = -^(2-w)rr"^(l) 



4 



^i")(«^,i) = 5r(«^,i) = o 

that gives, dividing the SR by the factor 2{w — 1), 
-L5{w) + {w + 1)Lq{w) 

= - E Wi'lliw) Cf ^"^(1) + % + 1)73 ri;)(^) rf (")(1) 
and using one gets finally : 



in). 



(71) 



(72) 



-L5{w) + {w + 1)Lq{w) 

= 2E ^1/1(1) - 4(^ + 1) E A^^;^ 4llM 41(1) . (73) 



6 Basic results. 

Let us recall the two sum rules that we have obtained in the two preceding 
sections : 

Uw) = -6 E ^E^yl A/lil) A%{w) (74) 
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-L5(w) + (w + l)Uw) 



w 



'3/2' 



w] 



(75) 



Due to the equations of motion, the functions Li{w) {i = 4,5,6) are not inde- 
pendent and, as shown in ref. 0, are given in terms of the elastic IW function C.i'w), 
a subleading function ^^{w) and the A parameter (A = — m^) : 



2 

w + 1 



(76) 



Therefore, from (fT^ - ffTUj) we obtain the interesting relations, valid for all w : 

n 



(77) 



[w 



+ l)EA^S^i/1(l)- 



'3/2' 



w) 



w] 



(7J 



These remarkably simple relations are the basic results of the present paper. 
They reduce to the known results (jl^ for t;; = 1. It is important to notice that 
both the subleading quantities AC,{w) and C,3{w) can be expressed in terms of leading 



quantities, namely the IW functions Tj^\w) and the level spacings AEj^' (j 

Very remarkably, equation ()77|) shows that the leading IW function ^{w) appears 
constrained to be a combination of the averages 



1 3^ 
2' 2' 



(j 



1 3 ' 
2 ' 2 / 



(79) 



or, conversely, the fundamental constant A is given by the ratio of functions 



A 



2{w + l)Y,^E. 



{") An), 



3/2 ^3/2(1) ^3/2 ('"^) 



+ 2Y.^E^%r^i;l{l)rt]l{w 



(80) 
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It is worth to underline that all other subleading IW functions except Q (1) 
cannot contribute to the l.h.s. of the SR (jSUj) because the polynomials (j^Tj) vanish 
at the frontier of the domain {wi,Wf,Wif) = {w,l,w). In the case of the B* SR, 
equation (jHOI), for the same reason only ^''''(1) and t^''\i) survive in the l.h.s. at 

{Wi,Wf,Wif) = (w, l,w). 



7 Phenomenological discussion. 

To illustrate our results for the subleading form factors, we will concentrate on 
some functions that play a role in the analysis of B ^ D{D*)£i', and about which we 
can get information from the relations obtained in the preceding Section on A^{w) 
and ^3{w). Of particular interest are the functions 

L,{w) and VM = Y^, (81) 

and their values and derivatives at zero recoil. 

The function L^lw) appears at first order in l/mq in the differential semi-leptonic 
rate of B ^ Diu fPj. This subleading IW function is specially important, but can 
be expressed, from (|7H| . in terms of C,{w), A and the commonly used function ri{w) 
(see for example ^2] and references therein) : 

Li{w) = -A^iw) [1 - 2ri{w)] (82) 

7.1 Check of the Bakamjian-Thomas quark models. 

We would like first to test whether the results found in the class of relativistic 
quark models of the Bakamjian-Thomas type are consistent with the sum rules 
found in this paper, in particular the w dependence. This class of models yield 
covariant form factors in the heavy quark limit that satisfy Isgur-Wise scaling, and 
Bjorken and Uraltsev SR. It is a class of models in the sense that one can choose the 
dynamics in the hadron rest frame, and then compute the corresponding Isgur-Wise 
functions with the boosted wave functions. 

The dynamics at rest that describes in the most accurate way the QQ, Qq and 
qq spectra (where Q and q denote respectively heavy and light quarks) is the phe- 
nomenological Hamiltonian set up by Godfrey and Isgur [22J, containing a confining 
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piece, a short distance piece with asymptotic freedom, plus spin-dependent interac- 
tions. 



7.1.1 Hypothesis of saturation by n = states. 

Using this model within the Bakamjian and Thomas scheme, one finds, for the 
IW functions of the n = states (n denoting the radial quantum number) 

aw) 



w + l 



W+l 

2 

w + l 



2 



A%iw)=T[%il)[-^) ^ (83) 



with 

p2 = 1.02 (84) 

for the slope of the elastic IW function, and, for the transition IW functions to the 
lowest P-wave states 

r^/°J(l) = 0.5394 42 = 1-50 

rJ/J(l) = 0.2248 ^2 = 0-83 ■ (85) 

It has been shown that these n = transition IW functions dominate the Bjorken 
[T], [2] and Uraltsev SR jS], that read, respectively 



p^ = ^+Ei-i^;^(i)p+2Eki/l(i)p 



E 141(1)1' -Ekl/1(l)P = I- (86) 

n n ^ 

Keeping only the n = states and the numbers quoted above we get 

l + |r(/°)(l)P + 2|r(;)(l)r = 0.882 

ki/1(l)P-|41(l)r = 0.240 (87) 

to be compared, respectively, with = 1.02 (j8l|l and with 1/4 in the r.h.s. of 
the second equation (|SB|). The n = states give a dominant contribution, and 
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saturate Bjorken SR at the 10 % level. Uraltsev SR is even more accurate in this 
approximation . 

Let us first test equation (jHUj) . saturating it with the n = states : 



A 



1 



Inserting the phenomenological IW functions ()83|) and the values for AEf'^ (j = 
|, |) obtained in the same BT scheme j23j : 

- ^^1/2 = 0-406 GeV (89) 

we get indeed a value of A that is quite stable in the whole physical region of 
B -> D*iiy, l<w< 1.5. We find 

A = 0.513 ±0.015 . (90) 

The stability of the result for A is quite remarkable, and results essentially from 
the function f{w) = {w + 1)t!^\,\{'^) the r.h.s., that has a slope /'(I) = very 
close to the slope of ^{w), ^'(1) = -1.02. 

On the other hand, one gets, from (fTTjl . (fTSjl . (jHTj) and the n = approximation : 

r/(l) = 0.380 r/'(l) = -0.006 ri"{l) = 0.0003 . (91) 

For the moment, let us notice that the accuracy of the n = approximation 
depends strongly on the considered quantity, at least at w = 1. We have seen that 
in Bjorken and Uraltsev SR the n = states dominate, but the precision is quite 
different in both cases. As for the subleading quantities, defining 



EA£;|;>|;^(i)l 

R = ^ (92) 
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one gets in BT models, keeping only n = 

Rbt = 0.174 . (93) 

7.1.2 Excited states (n ^ 0) contribution. 

Let us now discuss how this can be modified by ri 7^ states. Although the 
results obtained keeping only the n = states are encouraging, one must address the 
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question of the n 7^ excited states contribution to the SR. This question can have 
a clear cut answer within the BT scheme, but asks for further numerical calculations 
to compute Aii^]"'* and the ly-dependent form factors Tj^\w) {j = |, |), and will be 
done in a near future. 

At w = 1, we can for the moment state that a sum including higher n 7^ states 
leads, for the quantity (jHSI), to the value [TT] 

Rbt = 0.24 (94) 

giving 

r/(l) = 0.34 , (95) 

that differs from (j9T| by 10 %. Therefore, it is of importance to compute the 
contributions of n 7^ states for the different quantities, in particular to A of which 
()90|) is only a lower limit, due to the positivity of the different contributions. It is 
also worth to check whether the inclusion of the n 7^ states in ()80|) yields indeed 
a constant. 

In practice, one is anyway confronted to sums truncated to a definite Umax- It 
is of importance to notice that the dependence of the sum on rimax requires the 
consideration of the radiative corrections. On the one hand, one approach proposes 
to identify the renormalization scale fi with Ai?*^"""''"^) [22] • On the other hand, 
another point of view ^2], ^Bl distinguishes between the cut in n, given by a scale 
A such that Ai?^"""''"^) = A and the renormalization point /i, although eventually 
both scales can be chosen to be proportional. The discussion of the contribution of 
the higher n states is not simple and cannot be done without including the radiative 
corrections. 

7.2 Comparison with the QCD Sum Rules approach. 

Although more precise calculations remain to be done in the BT model, let us 
qualitatively compare with other approaches. 

The approach used up to now to obtain information on the subleading func- 
tions has been the QCD Sum Rules (QCDSR) approach |26. (for a review, see |EZ!)- 
Moreover, radiative corrections to the subleading l/rriQ corrections have also been 
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computed within this scheme in these works. For a recent discussion of the sublead- 
ing IW functions and their radiative corrections see and 

We must first notice that the resuhs obtained from the SR of the present paper 
and those of QCDSR are quite different in spirit. In our approach we have used the 
BT quark model to compute the r.h.s. of the SR (ffTjl and (|78|l . while in the QCDSR 
approach one computes directly the l.h.s. of these equations. 

The subleading IW functions are non-perturbative quantities, and their calcula- 
tion within the QCDSR approach is to some extend mo del- dependent because it is 
subject to a number of approximations. 

Hence the interest of having information on these non-perturbative quantities 
within the present method of Bjorken-like SR. We have considered here only a limited 
number of quantities, namely the subleading corrections of the current perturbation 
type : L4^{w), L^{w) and Lq{w) in the notation of Falk and Neubert p. Moreover, 
the radiative corrections to these quantities within the present approach have not 
been computed. We must compare the values obtained to the ones of the QCDSR 
without including radiative corrections. 

Without including QCD corrections, the QCDSR method gives the values [2Zj 

A = 0.50 

-,(1) = \ 

r/'(l) = (96) 

and sets [IS] 

^"(1) = (97) 

that are qualitatively consistent with our results (jHH)? (jSnj- 

Notice that, as already pointed out in [11^, the QCDSR algebraic value ri{l) = ^ 
would correspond to the value 

Rqcdsr = ^ (98) 

that is very close to the value including ri 7^ states [201, (HH]- 

As we realize in this comparison with the results of QCDSR, we have only com- 
puted in our approach a part of the subleading non-perturbative corrections, namely 
the 1/wiQ perturbations to the current. This is a part of a larger program that should 
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include the subleading quantities related to the perturbations of the Lagrangian, 
namely Li{w), L2{w) and L^{w) or, in the more usual phenomenological notation, 
Xi(tf), X2{w) and Xz{w) (see, for example [TH] and [2H])- 

8 Conclusion and outlook. 

In this paper we have shown that the consideration of the non-forward amplitude 
leads to powerful results for the subleading form factors at order I/tuq, at least in the 
case of the functions that correspond to perturbations of the heavy quark current, 
L4^{w), L^{w) and Lq{w), or equivalently A^(w) and ^3(1^), where Ow) is the leading 
elastic IW function and A = itlb — rrib. 

The parameters l/mb and l/rric are independent, and we have in this paper 
studied the sum rules coming from the terms in l/rrib, i.e. taking rric —>■ 00. The 
method in this case appears somewhat cumbersome but straightforward. The final 
results are very simple. We have considered in the SR intermediate states with 
orbital angular momenta of the light quark i = 0,1 and 2. The consequences that 
we draw from these states made explicit have been easily generalized to all £. 

Within the framework of the OPE and the non-forward amplitude, the sum 
rules of the type B{vi) D^'"'\v') B{vf) that depend on the three variables 
{wi,Wf,Wif) = (fj ■ v',Vf ■ v',Vi ■ Vf), allow to write A^{w) and (,3(10) in terms of 
leading quantities, namely the transition IW functions Tyl{w), t!^^1{w) and the 
corresponding level spacings AE["1 AE^yl This has been possible by taking the 
limit to the frontier of the domain {wi,Wf,Wif) = {w,l,w). Then, most of the 
excited intermediate states that contribute to the hadronic side of the SR vanish, 
because zero recoil is chosen on one side, namely Wf = 1, and only the P-wave IW 
functions Tyl{w), Tyliw) survive. As a result, the fundamental quantity of HQET 
A appears to be a ratio of leading functions and ^,3(10) is also given in terms of 
leading functions. 

To proceed further phenomenologically, we have used as an Ansatz for these 
functions the results of the Bakamjian-Thomas quark model, that gives covariant 
form factors in the heavy quark limit, satisfies IW scaling and also Bjorken and 
Uraltsev sum rules. One obtains in this way for a very wide range of w the expected 
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constancy for A, with a numerical value of the order of 0.5. This value is in agreement 
with the QCDSR approach. We find also numerical agreement between our approach 
for the ratio of functions ri{w) = C,3{w) / A^{w) and the QCDSR results. It must 
be emphasized that the confirmation of the results of the QCDSR approach by 
our rigorous method (with the phenomenological input of the BT model) is quite 
encouraging, since both methods are very different in spirit. 

The program to study the subleading form factors presented in this paper should 
be pursued in several directions. First, one should also compute the - in principle 
independent - sum rules of the l/rric type (c being the intermediate quark), that 
should be consistent with the ones of the l/m^ type computed here. We safely 
conjecture that no new SR for the form factors A^{w) and ^3(10) will be found. 
Second, one should compute within the BT scheme the n ^ contributions to the 
r.h.s. of the SR (|77jl and (f7H|) . Thirdly, one should study the other subleading 
quantities in the same spirit, namely the form factors that come from perturbations 
of the Lagrangian, i.e. Li{w), L2{w), L^{w), that remain rather uncertain in the 
QCDSR approach. However, it is not clear that usable relations could be obtained 
in this case, in terms of computed form factors in the BT class of relativistic quark 
models. But this direction should be pursued. Finally, one should follow and discuss 
the experimental situation for the n = P-wave IW functions rf'j^^iw), Ty\{w). As 
shown in Appendix B, the observed values of t[^1{1) do not seem to fit at present 
Uraltsev sum rule and the predictions of the BT quark models. However, the |^ 
states are very wide, and new semileptonic and non-leptonic data are necessary to 
finally settle the question of the values for Tyl{w). 
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Appendix A. 



We give here a proof of the identity between the subleading quantities S defined 
by Uraltsev P] and ^3(1) defined by Falk and Neubert 



^3(1) = S 



(A.i: 



Uraltsev expands the matrix element between two B* mesons for small velocity 
transfer u (formula (14) of the first ref. [H]), 



< B*{e',u)\QiDjQ{0)\B*{e,0) > 
-^u,ie'* ■e) + ^ [{e'* ■ u)s, - e';{e- u)] + 0{u 



(A.2) 



where Q is the heavy quark field. 

On the other hand, using the formula of Falk and Neubert (3.4) [9j, that is valid 
for all r° 

< B*{e' ,u')\QT''iD^Q{Q)\B*{e,u) >= -Tr \^o.{u,u')B*{e' ,u')r"B*{e,u)] (A.3) 



where 

one can write, in a covariant way 

< B*{e',u)\QtDjQ{0)\B*{e,v) >= -Tr [^j{u,v)B* {e ,u)B*{e,v) 
Computing in terms of the C,j{u,v) the matrix element (A.2) one has 

< B*{e' ,u)\QiDjQ{0)\B*{e,6) > = -Tr kj{u,v)B* {e' ,u)B*{e,v) 



with 



u 



[Vl + iP,u) v = {l,6). 



Therefore 



< B*{6',u)\QiDjQ{0)\B*(e,0) > = -Tr ij{u,v)f P+{u)P+{v 



—Tr 



{^+{w){u + v)j + ^^{w){u - v)^ - Uw)^j) tP+{u)P+{v)f (A.8) 



where w = y/1 + iP. 



(A.4) 



(A.5) 



(A.6) 



(A.7) 
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Taking into account that under the trace (A. 8) 



2e+(l)+e3(l) = 



(A.9) 



and the relation 



^-(1) - -I 



one obtains, after some algebra, expanding in powers of u, 



< B''{e',u)\QiDjQ{0)\B*{€,6) > 



(A.IO) 



A . . ^ ^3(1) 



[e'* ■ u)ej — e'*{e ■ u) 



that compared with (A. 2) demonstrates the identity (A.l). 



(A.ll) 



Appendix B. Comment on the experimental 
situation for ti/2(w), T3/2(w). 

The aim of the present paper has been a theoretical one. However, a brief com- 
ment on the experimental situation for the states D**{j^, J^) and the corresponding 
IW functions Ti/2{w), T-i/2{w) is in order. 

In the framework of the Bakamjian- Thomas quark model, we have given the 
prediction of the shape of the functions Ty\{w)^ ''"3/2 ('"^) n = states. We 

have seen that these n = IW functions almost saturate the Bjorken and Uraltsev 
SR, and give quite reasonable results for the subleading quantities A^(w) and iziw). 

However, for the |^ states, the present experimental situation seems at odds 
with these expectations. 

The more complete experimental analysis, from Belle, reports four excited states 
D** above the ground state in non-leptonic decays B~ D**^7t^ [31J, two broad 
and two narrow, that complete the expected number of n = states with parity 
P = +. The wide states should correspond to the (j^, J^) states Q^, 0+^ , Q^, 1+^ , 
and the narrow to (^|^, 1^^, Q^,2+^. In the following, we denote these states D-'j. 
Indeed, the strong decays proceed through Z^^^ Div, D^^ D*n (^-wave). 
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D-^' D*n, Dl' Dn, D*n (D-wave). The product of branching ratios 5(5" 
D**^n^) X B{D**^ D^*^^n^) have been measured. Assuming factorization of vr^ 
emission, assuming also that the states 1"*" are unmixed, and using a simple quark 
model for elementary pion emission for the decays B D**7r to estimate needed 
spin counting coefficients, one finds qualitatively, as we will see below, a magnitude 

\r[%{wo)\ - \4%{wo)\ (B.l) 

where wq = 7m^^^** is the value of w for = = 0. This situation is supported 
by older experiments measuring B D**7r ot B —>■ D**lv If the identification 
of these states is the correct one, the experimental situation for the states \^ is 
at odds with the expectation of the Bakamjian-Thomas quark model. However, 
an approximate saturation of Bjorken and Uraltsev SR by the = states is not 
excluded within la. 

However, before concluding about these states, new experimental confirmation 
is needed, mostly for the very wide |^ states. On the other hand, on the phe- 
nomenological side, one should take into account l/nic and l/m^ corrections in the 
decays B{B D**7r) 10^, and moreover a sensible theoretical scheme for the decays 
D** D^*'^T{ is needed, including l/rric corrections. New data on the semileptonic 
decays B D**ii>, that would allow to extract directly the functions Ty\{w), 
tI^^\{w) would also be very welcome. A needed detailed analysis to extract Ty\{w), 
tI^^\{w) from present data is beyond the scope of the present paper. Our aim below 
is to have only a qualitative estimation. 

The Belle data on the candidates to the four P-wave states D** are the following : 

D** (1^,2+) 

Ml'' = (2460.7 ± 2.1 ± 3.1) MeV 
tI'' = (46.4 ±4.4 ±3.1) MeV 
B{B- dI'' \-) X B{dI'' ° ^ D+TT-) = (3.5 ± 0.3 ± 0.5) x 10"^ 
B{B- dI'' \-) X B{dI'' ° ^ D*+n-) = (2.0 ± 0.3 ± 0.5) x 10"^ (B.2) 

D** (r,i+) 

Ml'' = (2423.9 ± 1.7 ± 0.2) MeV 
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rf ^ = (26.7 ± 3.1 ± 2.2) MeV 
B{B- dI^^ V) X B{D^^^^ ° ^ D*+7r-) = (6.2 ± 0.5 ± 1.1) x 10"^ (B.3) 

D** (1^,0+) 

Mq/^ = (2290 ± 22 ± 20) MeV 
rj/^ = (305 ± 30 ± 25) MeV 
B{B- ^ Dq/' %-) X B{dI^^ ° ^ D+TT") = (5.5 ± 0.5 ± 0.8) x 10"^ (B.4) 

D** (|+, 1+) 

M^^ = (2400 ± 30 ± 20) MeV 
r}/^ = (380 ± 100 ± 100) MeV 
B{B- dI'^ %-) X B{dI'^ ° ^ D*+7r-) = (4.1 ± 0.5 ± 0.8) x 10"^ . (B.5) 

S /2 

Assuming that these states decay essentially into two-body modes, i.e. B{D2 — > 
{D + D*)7r), B{dI'^ 0*71), B{Dy^ Dtt), B{D\'^ D*7r), the following 
branching ratios are given by a Clebsch-Gordan coefficient 

^p3/2 ^ ^*+^-^ _ 5(Dy' ° ^ D+TT") = B{dI^^ ° ^ D*+7r-) = - . (B.6) 

3 

To estimate 5(^2^^ ^ D+tt^) and B[Dy^ ^ — > D*+7r~), we use the spin 
counting of the non-relativistic quark model. 

In the limit of neglecting spin- dependent perturbations in the spectrum, i.e. as- 
suming the pairs {D,D*), (I?2^^, -D^^^) and (D}^^,!)^^) to be degenerate and that 
all the Dj are degenerate, simple angular momentum calculations give, for the total 
widths : 

r(D^/') = r(Df ') r(Dy') = T{Dy^) (b.7) 

r(Df ' ^ D*tt) - -r(D^/' ^ Dtt) . (B.8) 

This last relation gives the needed spin counting coefficient. 

It is easy to obtain this factor by realizing that to have the D wave (1 denoting 
the quark emitting a pion and taking Oz along the pion momentum) one needs 
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the operator ((T^/Ctt) 6xp(i2;iA;^) — > ikla^Zi. We have then, for the non- vanishing 
amphtudes 

M{D^^^ ^ Dtt) = < 1 0, 1 0|2 > < 0|a^|l > < 0|3^i^|l > 
^p3/2(±i) ^ ^*(±i)^) = < 1 0, 1 ± 1|2 ± 1 >< 1 ± l\al\l ±1x0 0|X|1 > 

(B.9) 

that gives 

^p3/2(±i) ^ ^*(±i)7r) = M(L>2/' ^ L>7r) (B.IO) 

and hence (B.8). 

We now take into account the actual masses. Since both D^^^ — > Dn and D2^^ 
D*7r proceed through the D-wme, we will have 

^^C^^--^-0.40 (B.ll) 

in an obvious notation. Therefore, we obtain the branching ratios 

5(£)3/2 _^ ^ 0.48 B{dI^^ ° ^ D*+n-) ^ 0.19 . (B.12) 

Prom these BR we find, adding the errors in quadrature 



B{B- dI^^ = (7.3 ± 1.2) X 10-^ (from L>f ^ ° ^ D+n') 
B{B- dI^^ \-) = (10.5 ± 3.1) X 10-^ (from ^ ° ^ D*+n-) 



(B.13) 



We realize that the value for B{B — > D^^^ °7r ) differs if one obtains it from 
£)3/2 ^ ^+^_ ^^^^ ^3/2 ^ 2?*+7r-, although they agree within la. Using 

(B.5) for the other modes, and taking into account the uncertainty from both results 
(B.13) one finds 

B{B~ L>f ^ %-) = (9.8 ± 3.8) X 10-^ 
B{B- dI''^ %-) = (9.3 ± 1.7) X 10"^ 
B{B- L>y^ %-) = (6.1 ± 1.4) X 10-^ 
B{B- dI'"^ %-) = (8.2 ± 1.4) X 10"^ (B.14) 
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The decays B~ — > Z}**°7r~ proceed through two different diagrams : a color- 
allowed diagram with vr" emission, and a color-suppressed diagram with _D**° emis- 
sion. We will now assume that the vr" emission diagram dominates, although this 
hypothesis could be incorrect, as we argue at the end of this appendix. However, 
assuming factorization of vr" emission and that the states 1"*" are unmixed, we find 
for the decay rates, from [21 : 

T. IT/ 12 r2 P 



V,,\' -J-\M{B ^ D**it)\^ . (B.15) 



167r 



\M{B ^ '7r)|2 = 2mD"mB(mB + mD**)\wl - lf\Ts/2{wo)\^ 



\M{B ^ = 2mD**mB{mB - m,,«)'(^o + 1)'K' " 1)^3/2(^0)!' 



\M{B ^ Dy\)\' = 4mD',mB{mB + m^«)'K - l)'|ri/2(w;o)|' (B.16) 

with 

^o = ^r — P=^ — B.17 

zmBTriB** 2mB 

the subindex denoting the value of w for = = 0, and m^.. the mass of the 
corresponding D{j^, J^) state. 

It is interesting to notice that the rates (B.15), (B.16) are given by the expressions 

2 



r(5 ^ DT^) = V{B Dr^) = f }^ 'l^y^^' 

T{B D\'\) = TiB Dl'\) = f^l^^f^^^ 



/l+r2 

(B.18 

(l+r 

(B.19) 

where r = '^^'^^i'^^ and r = respectively in the first and the second relations. 

TUB "triB ^ 

The equalities Y{B D^^'^n) = T{B Dl'^-n) and T{B Ol'^-n) = T{B 

1 /2 

Dq tt) follow from heavy quark symmetry, since there is a single helicity amplitude 
in all decays. 

Using the central values for the masses, but taking into account the errors in 
(B.14), we find respectively for the different modes, roughly : 
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B^Dl'^n |r3/2(l. 30)1 = 0.29 ±0.06 

fi^DpTT 1x3/2(1.32)1 = 0.27 ±0.04 

,2 (B.20) 
E^Di/V |ri/2(l. 33)1 = 0.35 ±0.04 

B ^ ny^n |ri/2(l. 37)1 = 0.39 ±0.06 

Within Icr there is consistency between the different determinations of |t3/2(wo)| 

and |ti/2(u'o)|, but errors increase considering both determinations. We conclude 

safely that we will have the numbers 

|r3/2(1.31)| = 0.28 ±0.06 

|ri/2(l. 35)1 = 0.38 ±0.07 (B.21) 

Extrapolating now with the Bakamjian-Thomas form factors (jEBI), ()85|1 for t^/2{w) 
and Ti/2(u^), we get 

1^3/2(1)1 =0.44 ±0.10 
|ri/2(l)| = 0.49 ±0.10 

to be compared with the values in the BT model |t3/2(1)|'^"^ = 0.54, |ti/2(1)|'^"^ = 
0.22. We find agreement for |r3/2(l)| within errors, but |ri/2(l)| is much too large 
compared with the BT model. 

Saturating Bjorken and Uraltsev SR with the n = states, we get a contribution 
to the Bjorken and Uraltsev SR that lies, within lo", in the following range, keeping 
only the n = states : 

p' = \ + ki/2(l)|' ± 2|r3/2(l)P = 0.90 ± 0.28 

|r3/2(l)|' - |ri/2(l)|' = -0.04 ± 0.18 . (B.23) 

Therefore, within lo", low values for are not excluded but the second value is too 
small compared to the r.h.s. ^ of Uraltsev SR (j86|l . 

We must keep in mind however that the estimation (B.22), that relies on the 
simple hypothesis of the dominance of vr" emission could be incorrect owing to two 
facts. First, in many decay modes the color-suppressed diagrams are empirically 
not so suppressed. Second, the diagram of D**^ emission is not computable on 
the ground of first principles in the BBNS QCD factorization scheme jHB], since 
the emitted meson is composed of heavy-light quarks. Therefore, one must keep 
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in mind that new data on semileptonic decays B D**iv (where the statistics 
is much smaller than in non-leptonic decays), that directly measure the functions 
Ti/2{w) and r3/2(u^), are necessary to settle the question of the magnitude of ri/2(l), 
r3/2(l) and their comparison with Bjorken and Uraltsev SR and with the predictions 
of Bakamjian-Thomas models. 
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